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We present a general method to derive the classical mechanics of a system of identical particles 
in a way that retains information about quantum statistics. The resulting statistical mechanics can 
be interpreted as a classical version of Haldane's exclusion statistics. 
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Particle statistics enters quantum physics in two re- 
lated but logically distinct way. The first one is related 
to the symmetry of the wave function, or more generally 
to phase factors associated with the exchange of identi- 
cal particles. The other is related to entropy, i.e. to the 
counting of quantum states, and is expressed through the 
Pauli exclusion principle and the phenomenon of Bose 
condensation. Some years ago, Haldane pointed out that 
it is possible to have a certain kind of quantum statis- 
tics of the second type, so called exclusion statistics, 
without any reference to wave functions or exchange fac- 
tors Q| . While the exchange phase factors characterizing 
fermions, bosons or (in 2 space dimensions) anyons [^]|| 
are intimately connected to quantum mechanics, there is 
no logical reason for not having effects of the second type 
of statistics even in classical systems. In fact, as is well 
known such effects must be put in by hand in order to 
avoid the Gibbs paradox in classical statistical mechan- 
ics. 

In this letter we show that it is possible to formulate 
a classical mechanics that builds in the effects of quan- 
tum statistics at the Lagrangian level. The dynamics, i.e. 
the equations of motion, will not depend on the statistics 
parameter, but the counting of states, and thus the sta- 
tistical mechanics, will. In the classical description the 
statistics is expressed through the occupation of phase 
space volume. Thus, each new particle introduced in the 
system will reduce the available phase space volume for 
the other particles, and the degree of reduction defines 
the classical statistics parameter. 

We first sketch a general formulation of the problem, 
and then present some specific results for two cases, non- 
interacting charged particles in a strong magnetic field 
and vortices in the Landau-Ginzburg-Chern-Simons the- 
ory. Both these systems are of interest for the quantum 
Hall effect. In both examples, we show that the result- 
ing classical statistical mechanics is a classical version 
of Haldane's exclusion statistics. Our examples are two- 
dimensional, but just as in the case of exclusion statis- 
tics, there is no reason in principle for our construction 
not to work in an arbitrary dimension. Below we will 
only present the general ideas and some of the main re- 
sults. A fuller account including calculational details can 
be found in ref . 01 . 



Consider a general quantum system and a subset of 
states \ip z ), which is labelled by a set of complex coordi- 
nates z = {z\, Z2, zn}. These may be the coordinates 
of a system of (identical) particles or the coordinates of an 
N soliton configuration. We only assume that the wave 
function evolves smoothly with a change of these coordi- 
nates, and that it is symmetric under an interchange of 
any pair of the N coordinates. To define the correspond- 
ing classical mechanics, consider the constrained system 
where the evolution of the full quantum system is pro- 
jected to the manifold, A4, defined by the (normalized) 
states |V>z)- The Schrodinger equation of the full system 
can be derived from the Lagrangian, 



L = ifi{ip\ip) - {ip\H\ 



(1) 



and the Lagrangian of the constrained system is obtained 
from this by restricting to the subset of states \tp z ). 
Expressed in terms of the coordinates z, it is of the 
generic form, 

L(z, z) = A-~Zi + A z Ji - V(zi,Zi) , (2) 

where A z is the Berry connection, 

A Zi =ih(^ z \dM , (3) 

V is the expectation value of the Hamiltonian in the state 

An important special case is when the state vectors 
which define M. are, up to normalization, analytic func- 
tions Of Zi, 



\ip z ) = AT(z,z)|(/) z 



(4) 



where \<j> z ) denotes the analytic part of the state vector, 
and J\f(z, z) is the normalization factor. The vector po- 
tentials are then given by 



ihd Zi lnAT(z, z) 



(5) 



and M. is a Kahler manifold where the Kahler potential 
is related in a simple way to the normalization factor, 



K(z,z) = h\n\Af(z,z)\- 



(6) 
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The geometry of M. is described by the field strength 

f ZtZj = dgtAj - d Zj A- = id- Zi d Z] K{z,z) , (7) 

in terms of which the sympiectic form, lu, and the met- 
ric, ds 2 , takes the forms u> — ~f ZiZj dz~i A dzj and ds 2 = 
—2if-. Zj dzidzj. The importance of the sympiectic form 
u> is that it determines the Poisson brackets, and thus, 
together with the Hamiltonian, defines the classical me- 
chanics H . The Poisson bracket is 

{^B} = f-i j d Zi Ad- z .B, (8) 

and the equation of motion can then be written as 



ii = {zi, V} 



(9) 



To be more specific we now consider a system of 
charged particles moving in two dimensions in the pres- 
ence of a strong magnetic field that restricts the avail- 
able states to the lowest Landau level. In this example 
we can explicitly derive the metric and sympiectic form, 
and show that they can be obtained from a Kahler po- 
tential. For calculations, it is convenient to consider par- 
ticles moving on a sphere. On a unit sphere penetrated 
by 2j units of magnetic flux a particle with unit charge 
has a total angular momentum J = j + L, where L is 
the orbital angular momentum, and the lowest Landau 
level corresponds to L = with a 2j + 1 degeneracy [|| . 
We shall use the notation of ref. [Q and define a coherent 
state by rotations of a minimum uncertainty reference 
state |0), 



= D(z)\0)=e 



*J+ e vJo e - 



10) 



(10) 



Here the complex coordinate z is defined via a stere- 
ographic projection, and the rotation operators, D(z), 
form a unitary and irreducible representation of the ro- 
tation group, generated by J m , rj = ln(l +zz), and |0) is 
annihilated by J+. 

Fully symmetrized and antisymmetrized states corre- 
sponding to fermions and bosons are given by 



|z,±) =W(z,z) 



1 



r]$e Zi r J +\0) 



m V 



(11) 



with rj p the appropriate sign for the permutation P. The 
normalizations of these states are readily obtained from 
the properties of the D(z):s, 



iAA( Z ,z)r 2 = 5> P ;Q(i+^) 
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(12) 



For the case of N coinciding bosons, Zi — z, we immedi- 
ately get the following Kahler potential 



K(z,z) = Nh2jln(l + zz) , 



(13) 



and the corresponding metric, ds 2 = 2Nh/(l + zz) 2 dzdz, 
is just N times that of a sphere. 



To assess the effect of statistics in the classical descrip- 
tion we calculate the TV-particle phase space volume. Fol- 
lowing Manton M, and Samols we can use ( |l3| ) to 
obtain the iV-particle volume from the volume of N co- 
inciding bosons, and the result is, 



V B 



1 



(A) 



N 



with A as the volume of the single-particle space, 

a t f ,_. MirR 2 

A = h u = h2] = —— = e$ , 

J sph ' 



(14) 



(15) 



which is h times the number of flux quanta (f>o = h/e that 
penetrate the sphere. Thus, for bosons the only effect of 
the indistinguishability of the particles is the factor l/Nl, 
and there is no further reduction in phase space volume. 

One should note that the classical phase space defined 
as above is everywhere a smooth manifold. This is dif- 
ferent from the case when the N-particle space is defined 
as a product of single particle spaces with identification 
of equivalent configurations. In the latter case the points 
corresponding to a coincidence of two particle positions 
are singular. Also note that the factor l/Nl here appears 
naturally from the geometry, not through any additional 
assumption about identification of points. 

For fermions a similar but technically more involved 
calculation can be done. The resulting phase space vol- 
ume is, 



V F = j ri (A-(N-l)hf 



(16) 



Compared with the Bose case there is an additional re- 
duction of phase space. The available phase space for any 
particular fermion is reduced with an amount h by each 
of the other particles present in the system. This can be 
understood as a classical version of the Pauli exclusion 
principle, and is consistent with the usual semi-classical 
interpretation of quantum mechanics, where each quan- 
tum state is associated with a phase space volume h. 
Note that there is a maximum number of particles al- 
lowed, N = 2j + 1, in which case the phase space volume 
( |l6| ) vanishes. This corresponds to all states with zero 
orbital angular momentum being filled, i.e. to a filled 
lowest Landau level, which is an incompressible state. 

The calculation can be done also in the case of anyons 
although there is no explicit expression for the normal- 
ization constant corresponding to (12). For details we 



again refer to [Q and only quote the result 



v(N-l)h) 



N 



(17) 



where the exchange phase of the anyons is uw. 

The expressions we have found above for the iV-particle 
phase space volume demonstrates a classical fractional 
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exclusion principle. Thus, each new particle introduced 
in the system will find the available volume reduced by 
a = vh relative to the previous one. The quantity a, 
i.e. the reduction in phase space volume, can be taken 
as defining the classical statistics parameter of the parti- 
cles. In the present case it is simply the (dimensionless) 
quantum statistics parameter v multiplied with Planck's 
constant h. 

We next consider a classical field theory with soliton 
solutions, namely the Chern-Simons Ginzburg-Landau 
(CSGL) theory, originally introduced as a field theory 
for the quantum Hall effect JlO| . 

L = J d 2 x[iH*D <p - ^im 2 - ^(M 2 - Po f (18) 
+ fihe^ p a fJ _d u a p \ , 

where is a complex matter field, a M a Chern-Simons 
field, m a mass parameter, A the interaction strength, 
Po the preferred density of the system and p a statis- 
tics parameter. For the original Laughlin quantum Hall 
states described by the model the statistics parameter 
takes the values p = l/[47r(2fc + 1)]. This theory has vor- 
tices (quasi-particles) as soliton solutions. In a certain 
approximation the dynamics can be described in terms 
of vortex coordinates alone, and a phase space description 
can be derived from the full theory. Again it is possible 
to calculate the phase space volume corresponding to a 
iV-vortex solution. 

The vortex configurations can be parameterized by a 
set of vortex coordinates z = {zi, Z2, z n }, just as the 
charged particles in a magnetic field discussed above. 
The precise form of the multi-vortex configurations for 
given coordinates is not known, but for critical coupling 
(A = 1) the existence of 7V-vortex configurations with ar- 
bitrary positions can be deduced With the system 
constrained to the manifold of A^-vortex configurations, 
a classical mechanics follows with a kinetic term for the 
iV-vortex system corresponding to a phase space with 
Kahler metric jjj. By use of earlier results obtained by 
Manton for the related relativistic abelian Higgs model 
H we find for the phase space volume of N vortices, 

V N = ^(A-4Trph(N-l)) N , (19) 

where the classical statistics parameter, as determined by 
the reduction in available phase space due to the presence 
of other vortices, is a = Airph = gh. We can interpret g, 
the classical parameter divided by h, as the dimensionless 
quantum statistics parameter. The value g — iirp agrees 
with the value of the statistics parameter as determined 
from Berry phase calculations with Laughlin wave func- 
tions p2[, or from the properties of vortices in the CSGL 
model |ic|] . 

We now discuss the statistical mechanics of the clas- 
sical systems just derived. Both these systems have the 



special property that the energy does not depend on the 
state, but only on the number of particles. This moans 
that the statistical mechanics is determined by the phase 
space volume Vjv, which has been determined in the pre- 
vious sections, and by the energy En. The classical parti- 
tion function is simply the total number of states, Vjy /h N 
multiplied with the Boltzmann factor, i.e. 

Z N = ?$e-0»" , (20) 

Using standard thermodynamics we get for the entropy 

S = N\a(l — ap) + NItl — — NItlN + N, (21) 

where a — vh or gh, and where we have introduced the 
classical phase space density p — N/A and neglected the 
difference between N and N — 1, which is irrelevant in 
the thermodynamic limit. 

In the systems we have considered the real two- 
dimensional space where the particles or vortices move 
is proportional to the phase space. Defining the pres- 
sure as, P = — (8F/dA) T , where A = V\ is the phase 
space volume for a single particle, we then get, by use of 
standard thermodynamic relations, the equation of state 

pP = p/(l-ap). (22) 

This expression shows that there is a maximum density 
p = 1/a allowed by the system, which corresponds to 
an infinite pressure and therefore to an incompressible 
state. For the phase space volume this means Vn = 0, 
i.e. there is no available phase space volume for any new 
particle added to the system. For the anyon system this 
situation corresponds to a completely filled Landau level. 
What is unusual about this is that the blocking, which 
can be interpreted as representing a generalized Pauli 
principle, shows up not only in the quantum, but also in 
the classical description of the system. 

Finally we show that the thermodynamics just derived 
can be viewed as a classical limit of Haldane exclusion 
statistics 

The statistical mechanics of particles with exclusion 
statistics can be derived by assuming that the total en- 
ergy can be written as a sum of single-particle energies 
and using the prescriptions for statistical weight given by 
Haldane separately for each single-particle energy level 
|l3| |l5| . The result for the entropy is 

S = + (1 - g)n k ] ln[l + (1 - g)n k ] (23) 

fc 

+ (1 - gn k )ln(l - gn k ) - n k lnn k } , 

where the sum runs over single-particle energy states, 
and g is the exclusion statistics parameter. D k is the de- 
generacy of the A:-th level and the quantum distribution 
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function is the average occupation number of the state 
k. 

Since each quantum state occupies the phase space vol- 
ume hP , with 22? the dimension of the single-particle 
phase space, we can relate n and p in the semiclassical 
limit by n — phP . In the Boltzmann limit, h — > and 
n — ► 0, all dependence on g in (^3|) goes away. If we, 
however, define the classical physics by the limit h — ► 0, 
g — ► oo and g/i 13 — ► a, where a is interpreted as a clas- 
sical statistics parameter ( p^ ) gets the nontrivial limit 
of 

5 = ^ D kh v [pk ln(l - ap k ) - pk hi(p fe /i) + p*.] . (24) 

If we further specialize to the case of fully degenerate 
states in a two-dimensional phase space, where the sum 
is simply replaced by the total number of available single- 
particle states, G = A/h, and where pk is replaced by 
N/A, we exactly regain (pl|). This demonstrates that 
the classical statistical mechanics discussed in the previ- 
ous section can be regarded as a special limit of exclusion 
statistics, different from the Boltzmann limit. Starting 
from the the equation of state for exclusion statistics par- 
ticles with the same energy, we can also derive ( p2| ) if we 
identify the the physical volume V with the one particle 
phase space volume A. 

One can also obtain exact results when the particles 
move in an external harmonic potential. Again one 
finds equivalence between the statistical mechanics de- 
rived from the classical mechanics of identical particles 
with a statistics parameter and the statistical mechan- 
ics derived from exclusion statistics in the classical limit 
discussed above Q]. 

Although the phase space in the examples discussed 
in this paper are two-dimensional, there is no obvious 
reason for the construction not to work in higher dimen- 
sions. One point of particular interest to study further 
is how the "classical fermion" theory discussed here can 
be applied to systems of interacting fermions moving in 
two or three dimensions. 
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